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Abstract 

A new statistical model for multiparticle production in e^e^ annihilation is 
proposed based on the idea of the longitudinal phase space with limited transverse 
momentum. The longitudinal rapidity space is divided into cells of equal size in order 
to take into account the Bose-Einstein correlations (BEC) with a finite correlation 
length 6y. The maximum entropy method is used to determine the probability 
distributions of the final state pions (or p mesons) for a given mean multiplicity, 
mean transverse momentum and mean total energy. Event simulation based on our 
model is performed in two extreme cases, a vr-model and a p-model. The vr-model 
assumes that only pions are produced directly. On the other hand, the /9-model 
assumes that only p mesons are produced directly and they decay into pions. A good 
overall fit to experimental data from PETRA to LEP energy regions is obtained for 
6y ranging from 0.6 to 1.2 in the vr-model. We found that the BEC plays a very 
important role to reproduce various correlation data. In some correlation data, 
resonance effect and conservation laws are also important. 

PACS: 12.40.Ee, 13.65.-Fi, 02.70.Lq, 05.40.-Fj, 25.75.Gz. 

1 Introduction 

With increase in energies of accelerators, accurate data on high-multiplicity events at high 
energies have recently been provided. This makes it possible to analyze correlations and 
fluctuations such as the Bose-Einstein correlations (BEC) |]T| and the intermittency 0, |^ 
in detail. These correlations and fluctuations of produced particles are rather new tools 
to study multiparticle production. By analyzing them, one may be able to extract useful 
information on the space-time size of the particle production region and the production 
mechanism |jl|]. Correlations and fluctuations may also provide a clear signal of quark gluon 
plasma formation expected in high energy nucleus-nucleus collisions. On the other hand, 
multiparticle production is a phenomenon observed commonly in various reactions such as 
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e'^e~ annihilations, lepton-hadron deep inelastic interactions, hadron-hadron interactions, 
and nucleus-nucleus collisions. A large amount of correlation data has been accumulated 
on various reactions. It has been found that there are some universal characteristics 
which are independent of the reaction type while some other properties depend on it. 
Phenomenological models which can describe some crucial properties of many particle 
correlations will be very useful for systematic investigations of those characteristics caused 
by correlations and fluctuations. 

The BEC have lately attracted considerable attention, in particular, because they may 
give information on the space time structure of hadronic source and they may cause large 
fluctuations of particle density in the phase space. Although many existing phenomeno- 
logical models have not taken into account the BEC completely yet, there are some models 
(or event generators) which allow computing the two identical particle correlation func- 
tions. In most cases, the BEC are generated more or less artificially. For example, they 
are calculated by using (i) the Fourier transformation of the model's source function|^ 
or (ii) the Wigner function]^, (iii) by weighting every event or (iv) by modifying the 
distribution of momentum difference of identical particle pairs in each event |^. However, 
some of those models cannot describe the BEC in a single event. Motivated by those 
consideration, we would like to propose a new statistical model. It is constructed by the 
maximum entropy method |^ and two versions of the model, the vr-model and the 
p-model are applied to e^e^ annihilation in this paper. The most characteristic point 
of our model is that the BEC is taken into account in quantum statistical level with a 
characteristic correlation length {5y) defined in the rapidity space. 

This paper is organized as follows. In Sec. I, our model is explained and the distribution 
functions of the final state hadrons are derived. Sec. 11 is devoted to the explanation of 
the simulation method including the procedure of parameter determination. Comparison 
of results obtained from simulation with the experimental data from PETRA to LEP 
energies is given in Sec. IV. Concluding remarks and discussions are given in Sec.V. 



2 Statistical model based on Maximum Entropy Method 

We first consider the vr-model, i.e., we assume that only pions are produced according to 
a statistical distribution. A system consisting of many pions produced in a single event 
may be decomposed into three subsystems, which consist of like sign pions, that is, vr"*", 
TT" and tt" subsystems. For example, consider a vr"*" subsystem in the longitudinal phase 
space with limited transverse momentum[0]. The longitudinal axis may be identified 
with the direction of the initial quark and antiquark produced by the virtual photon in 
center of mass system. Then the rapidity space in the longitudinal phase space is divided 



into many cells of equal size 6y^ and we consider a probability to find tt"'" mesons in 
each cell. The probability of finding 7r"'"s in the i-th cell, P/"'"* is normalized as 

oo 

E = 1- (1) 

ni=0 

The summation index rii runs from to infinity because of the Bose-Einstein statistics. 
The mean vr"*" multiplicity and the mean energy of the tt"*" subsystem are given, respectively, 
by: 

K) = EE^^^^^ (2) 
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and 

{E^)=J:J:^^^^Pt'^■ (3) 

i rii 

Here is the energy of a vr"*" in the i-th ceU: 

e . = mT cosh yi , (4) 
l/i = Vmin + {i- ^)Sy , (5) 



ymin 111 



(6) 



where i/i is the central rapidity value of the z-th cell, ymin is the kinematical minimum value 



of rapidity, = y {Pt^ + m-l is the mean transverse mass, and ^/s is the total center of 
mass energy. The same argument applies to both tt" and ir^ subsystems. Therefore the 
mean multiplicities and the mean energies of those subsystems are given as 

= (no) = = (n,,)/2, (7) 

and 

{E+) = {Eo) = = v^/3, (8) 

where {rich) is the mean charged multiplicity. The most probable distribution p/"''' for the 
TT"*" subsystem can be determined by the maximum entropy method which is well-known 
in information theorypj. According to a statistical picture, particles will distribute in the 
phase space in such a way that every possible state is realized with the same probability. 
This method has been applied in many fields in science including multiparticle production 
phenomenology. First, we define the missing information entropy 

5 = -EE^/"^i^^/"^- (9) 

i rii 

This S has the maximum value when p/"'"* is a constant, so-called equal a priori prob- 
abilities, which does not depend on i and rii. It is noted here that the probability P/"''' 
does not depend on pj"^^ when (i 7^ j). Therefore a probability of finding rii 7i"^s in the 
i-th cell and rij vr+s in the j-th cell simultaneously, p/"""'^'*, is given by 



The entropy S should be maximized with the constraints which correspond to the 
normalization eq.(|Tp, the mean multiplicity eq.(^ and the mean energy eq.(^) of the 
subsystems. For that purpose we consider F defined below instead of S by introducing 
Lagrange multiphers X\, A2 and A3; 

F= - EE^/"'^lnP/"') 

i rii 

i rii 

+ A2(EE^^^^^-K)) 

i rii 

+ HT.T.n.e.Pt'^ - {E^)) . (11) 
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By requiring that the variation of F is vanishing, one obtains 

P/-) = ^- . (12) 

Here, the Lagrange multiphers A*]^, A2 and A3 are rewritten in terms of the partition 
function Zj, the "partition" temperature T and the "chemical" potential [i: 

= 1 - e^^-^)l^ = i/Zi , (13) 

A2 = /i/T , (14) 

A3 = -l/T. (15) 

The parameters T and yU are determined uniquely for given and and they play 
a dominant role to reproduce the gross features of data on single particle spectra. The 
probability distribution functions for vr" and vr" subsystems are the same as P/"') for 
TT"*" subsystem. The parameter 5y is determined by fitting to data on the second order 
fluctuations of multiplicity. 

In order to study the resonance effect, we also consider another extreme case called 
p-model where it is assumed that only p mesons are produced in the same way as pions in 
the TT-model. In the p-model, one must take the spin degree of freedom into account. A p 
meson system is decomposed into nine subsystems: p+t, p^— p"*"!, p°t; P°~^) P°i) p T; 
p~— and p~|, where j, — and [ denote the three eigenstates of the spin = +1, and 
-1, respectively. The partition temperature Tp, the "chemical potential" pp and the cell 
size 6yp in the p-model are determined so as to reproduce the mean charged multiplicity, 
the total energy and the second order fluctuations of multiphcity. (The decay process of 
p meson will be discussed in Sec. 11 B) 

3 Method of simulation 

The method of generating events based on the 7r-model and the p-model is explained in 
this section. 

3.1 TT-model 

Events generating procedures are divided into three steps. The first step is determination 
of the parameters p and T. The second step is 6y search with trial events. By using these 
determined parameters, the final simulations are executed in the third step. 

Step 1: Determination of the parameters p and T 

To determine the values of T and p, we solve the following simultaneous equations for p 
and T with a fixed trial value of 6y: 

gni(/i-ei)/T 

(nc/,)/2 = 7 

= E i_g(,-..)/T > (16) 
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^ni{ix-ei)/T 



g.e(M-e<)/T 

^ 1 _ e(M-ei)/T' 

where observed value is used for the mean charged multiphcity {rich)- It should be noted 
here that fi and T are determined as functions of i.e. ii = iJ>{Sy), T = T{dy) . 



Step 2: Generation of trial events and determination of the pEirameter Sy 

1. In order to avoid a possible artifact due to a particular cell location, we now take the 
central value yi of the i-th cell as yi = ymin + — \ +f)^y-i where r is a homogeneous 
random number between and 1. 

2. The TT"*", 7r° and tt" subsystems are generated in all cells in the rapidity space 
according to the probability distribution P^^"'^ with a trial 6y^ l~^{Sy) and T(Sy). 
When a pion is produced in a certain cell, its rapidity is assigned by using a random 
number so that pions distribute homogeneously in that cell. This smearing of pion 
rapidity causes slight changes of {rich) and ^/s in comparison with the result of 
the step 1. We readjust the values of and T to reproduce correctly the mean 
multiplicity and the "mean" energy of the system. 

3. Charge conservation is demanded event by event. Events are discarded if the charge 
conservation is violated. 

4. The transverse momenta are generated according to the following distribution: 

^ (X e-^*'-/^-^ . (18) 
apj. 

Here (pr) denotes the observed mean transverse momentum at each y/s. 

The azimuthal angle (p of the transverse momentum is distributed at random be- 
tween and 27T. The transverse momenta generated in this way are assigned to all 
pions independent of their rapidities. A four momentum of a pion is then given as 

p'^ — ( ml. + cosh y, p^cosf/?, 

PTsimp, ^ ml + p^ sinh y ). (19) 

5. Approximate energy-momentum conservation is also required. For energy conser- 
vation, we set an energy window with width ±5E around ^/s. For momentum 
conservation, we also set the longitudinal and the transverse momentum windows 
with widths i^P^and 6Pt, respectively, around total momentum. The total en- 
ergy Etot, the longitudinal component of the total momentum Pj, and the transverse 
component Pt are evaluated in every event. An event is adopted only when it sat- 
isfies the inequalities 

^-SE< Et,t=Y.e <^ + SE, (20) 

all TT 
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and 

-6Pl< Pl=J2Pl < ^Pl, (21) 

all TT 

P| =1 E Pt \'< . (22) 

all 7r 

As an example, the distribution of Etot and the apphed energy window for -^5=34.5 GeV 
are shown in Fig.l. One has to take the window size as small as possible to keep 
the quality of conservation law while one would like to take it as large as possible 
to save the computation time. Thus one has to determine an optimum value of 
window size. We decided to take 5E = 0.2y/s and (5Pl = 0.2praax, where Pmax is 
the maximum value of the cm. momentum carried by a pion. By the way, we 
found that some physical quantities, e.g., the dispersion of the multiplicity distri- 
butions are rather sensitive to 6E when SE is large. We have confirmed that this 
strong (5i?-dependence disappears when SE = 0.2y/s, which implies that the physi- 
cal observables are insensitive to 6E if 6E < 0.2^/s. For p^-conservation, we use a 
multiplicity dependent window size 6Pt = \/n^ + uq -\- (p^) . 

6. The jet axis, thrust or sphericity, is determined event by event. New four compo- 
nents of a momentum vector of a pion are calculated by referring to this jet axis. 

7. The value of the parameter 6y is determined by fitting to appropriate experimental 
data which are sensitive to 6y. We have used the rapidity interval dependence of 
DlhK'^chf for the fitting. Here, D^h = \/(«c/x) - («c/i)^- 



Figure 1. 



Step 3: Final event generation 

A sufficiently large number of events are generated using the probability distribution p/""'' 
with /i, T and 6y determined in the preceding steps. 

Final results of those parameters are shown in Table I and the results of the fitting 
are shown in Fig.2(a)-(c), where ycut is the half size of the rapidity interval [—ycut,ycut]- 



Experimental data are taken from refs.|13, 14|. As shown in Table I, 0.6, 0.8 and 1.2 



are chosen as the values of the 6y at = 14.0, 34.5 and 91.2 GeV, respectively. In the 
determination of the best fit value 6y, we have put a special weight on the data point 
for ycut = ymax in order to reproduce the multiplicity distribution for full phase space as 
precisely as possible. It should be noted that ycut dependence of the D'^,^/{nch)'^, a kind 
of 2nd order fiuctuations of the particle density in the rapidity space, is reproduced well. 



Figure 2. (a)-(c) 
Table I 
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3.2 y9-model. 



The parameter Tp, fip and Syp in the p-model are determined in a way similar to the tt- 
model. The results are given in Table E. 



Table H 



Now we discuss the decay of p mesons e.g. — > n'^n^. Rapidities have to be assigned 
to each pion produced from one p meson of some rapidity Up- In the rest frame of the 
p meson, the maximum and the minimum rapidities of the two pions produced by the p 
meson are given by ±yd where 



Ed^^pl + ml , pd = ^M^/A 



ml 



and Mp and m^^ are the masses of the p and tt mesons. We use the following Breit-Wiger 
form f{Mp) as the invariant mass distribution of p meson: 

where mp=770 MeV and r=150 MeV. 

To impose the Bose-Einstein correlations between like charge pions, we use an "en- 
hancement factor" method. When a parent p meson has a rapidity y^, we divide the 
longitudinal rapidity interval [yp — yd,yp + ya] into K cells of equal size Sy^^. Here we take 
Sy^, = 6y. In order to determine the cell to which a pion, say tt"*", belongs, we introduce 
the enhancement factor: 

for i-th cell, (25) 



K + N+ 



where nf and A'"+ are the number of 7r'*'s already produced in the i-th cell and that in 
the interval [yp — yd,yp + yd], respectively. Transverse momenta pj,^ and pj,^ are assigned 
to the two pions, in the same way as in the tt- model. To reproduce the invariant mass 
of the parent p meson Mp, the azimuthal angle difference, a0 = 0i — 02, needs to satisfy 
the following equation: 

Mp = 2ml + 2 (6162 -PliPl2 - \Pti\\Pt2\ cosa0), (26) 
= + coshy^, (27) 
PLi = \/ml + \p^.\ sinhyi for i = 1 or 2. (28) 



4 Comparison of results with experimental data 

In this chapter we compare the results of event simulation with experimental data. We 
investigate effects of charge and energy-momentum conservations and also the resonance 
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effect by comparing the results from vr- and p- models. 20000 events of multiparticle 
production in e^e~ anniliilation liave been generated for both the vr-model and the p- 
model at each energy ^/s = 14.0,34.5 and 91.2 GeV. We use the experimental data 
reported by the TASSO Collaboration IH, |15|, |T|, [l^ for = 14.0, 34.5 GeV, and by 



the DELPHI |4|, |T8|, [T9|, |2g and ALEPH Collaborations 121|, ^ for 91.2 GeV. We present 
all results obtained from the vr-model. When significant differences are found between the 
two models, we present also the result of the p- model. 



4.1 Multiplicity distributions 

Charged multiplicity distributions P{nch) are shown in Figs. 3 (a)-(c). The experimental 
data||13|, O] are well reproduced by our vr-model. This is reasonable because the mean 



multiplicity and dispersion of the multiplicity distribution are used in the determination 
of the model parameters. It should be noted, however, that our model is able to predict 
the detailed shape of the multiplicity distribution at each energy. We also investigate the 
multiplicity distribution in terms of KNO variable. Our model gives good scaling behav- 
ior for both the full phase space(see Fig.4(a)) and small rapidity windows(not shown). 
The i/cnrdependence of the multiplicity distribution is reproduced correctly as shown in 
Fig. 4(b). Note that our results are approximately consistent with the negative binomial 
distributions. 



Figure 3. (a)-(c) 
Figure 4. (a) and (b) 



Dependence of multiplicity distributions on Sy has already been demonstrated in 
Figs. 2 (a)- (c). Decreasing 6y reduces the dispersion of the multiplicity distribution. Con- 
servation laws also affect the multiplicity distribution. In addition to "full conservation" 
event for which both energy-momentum and charge conservation are imposed, we have also 
generated events with "only charge conservation" , "only energy-momentum conservation" 
or "no conservation" by imposing only charge conservation or only energy-momentum con- 
servation or none of them in order to investigate the effects of conservation law^. It is 
found that the charge conservation causes the broadening of the multiplicity distribu- 
tion, while the effect of energy conservation is opposite as shown in Fig.5. The effect of 
energy-momentum conservation is stronger than that of charge conservation. 

We have investigated the charge correlations due to charge conservation. For this 
purpose, we consider the rapidity window dependence of the quantity Dlf^/2D\, where 

D"]^ is the dispersion of the vr"*" multiplicity distribution, ^J {n\) — {n+)^. The quantity 
D'^f^/2D^ can be written as 

Dch ^ ((n+ + n-)^)-(n+ + n-)^ 
2Dl 2Dl 

^ Dl 

^ It should be noted here that both energy-momentum and charge conservations hold on the average 
even in events with no conservation. 
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where use was made of the equahty (n^) = If there is the maximum correlation 

= n_ due to charge conservation, one has Dlf^/2D\ = 2 because {n^nJ) = {n\). On 
the other hand, if there is no correlation between and {n^n_) is reduced to (?^+)^ 
and hence Dl^/2D']^ = 1. Therefore, it is obvious that D^f^/2D'^ = 2 for i/cut = Umax- On 
the other hand, one can expect that Dlf^/2D']^ = 1 for very small i/cut because the corre- 
lation due to charge conservation will be maximally weakened when i/cut is small. Those 
expectations are indeed realized by our simulation as shown in Fig. 6. The correlation 
measure Dly^/2D']^ increases linearly from about unity as Ucut increases and then tends 
to saturate at the value 2 for ycut >3.5. The slope of the hnear rise may be regarded 
as a measure of local charge conservation. Unfortunately, it appears that corresponding 
experimental data are not available at present. We therefore urge experimentalists to 
provide such data. 



Figure 5 
Figure 6 (a) and (b) 



4.2 Rapidity distribution 

The single particle rapidity distributions are calculated and the results are shown in Fig. 7. 
Our model reproduces the experimental data[|I^ except for the observed dip structure at 
?/ = 00. In particular, the energy dependence of the central height is reproduced correctly 
by our model. Although our model produces a dip at |/=0 in qualitative agreement with 
experimental data, it is too deep and too narrow. Furthermore it tends to diminish as 
^/s increases. This tendency appears opposite to what is observed experimentally. This 
situation does not change even if one uses the p- model. The discrepancy may be due 
to the fact that three jet events are not included in our mo deig. By the way, we found 
that the dip is not reproduced if one uses the sphericity instead of the thrust to define 
the jet axis. We have also confirmed that the rapidity distribution is insensitive to both 
conservation laws and the value of 5y. 



Figure 7. 



4.3 Rapidity dependence of the mean transverse momentum 

There are interesting experimental data on the rapidity dependence of the mean transverse 
momentum, (pr(y)) HH- The observed (pt(|/)) is almost constant in the central region 



while it decreases as y approaches the kinematical limit as shown in Fig. 8 (a). On the other 
hand, transverse momentum of a pion is generated according to the exponential formula 
(18) with ^/-independent mean transverse momentum (p^) in our model. Therefore, it 
is very interesting to see if our model can reproduce the observed y-dependence of the 

^ The result of our model is systematically larger than TASSO data at -^=14.0 and 34.5 GeV, because 
integration of the dn/dy data reported by TASSO is smaller than the reported mean multiplicity. 

This dip structure had been investigated by TASSO collaboration[^. They report their analyses 
using the Lund-string model[E3| and the independent jet model |24[. 
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mean transverse momentum. The result of our model is shown and compared with data 
in Fig. 8 (a). Our model reproduces correctly the observed y-dependence of the mean 
transverse momentum. Therefore, one can conclude that the observed y-dependence is 
totally due to kinematical reasons. To confirm this conclusion, we have calculated {priy)) 
for events generated without energy-momentum conservation and/or without assigning 
the thrust axis. The result for 34.5 GeV is shown in Fig. 8(b). It is clearly seen that 
the energy-momentum conservation is most responsible for reproducing the correct y- 
dependence but the use of thrust axis has also some effect. 



Figure 8. (a) and (b) 



4.4 Bin size dependence of the scaled factorial moments 

We study the scaled factorial moments in order to investigate the fluctuations of particle 
density in the rapidity space. The following two kinds of scaled factorial moment are used 
for comparison with experimental data. 

The factorial moments are evaluated for the multiplicity of particles produced in the 
rapidity window yo < y < yo + AY. The rapidity window is divided into M bins of equal 
size AY/M. When N particles are produced in the window while particles are found 
in the m-th bin, the 'exclusive' scaled factorial moment is defined as 

=^ \h N{N-l)...{N-^ + l) /' 

where ( ) means an event average with fixed N. The 'inclusive' scaled factorial moment 
used in many literatures 0] is given by 

Fr = M-V E ^"'^^"^ ~ ~'^^^ ), (30) 

\ m=l ) ' 

where ( ) now means an event average for distributed A^. Both and F/" are equal to 
unity if there are only statistical fiuctuations( a binomial distribution for Ff^ or a Poisson 
distribution for -F/"). On the other hand, the Bose- Einstein distribution (eq.(12)) gives 
a typical example of non-statistical fluctuations. It is the multiplicity distribution in a 
single cell of size 6y in our model, corresponding to M AY/6y. In this case the inclusive 
factorial moments of the second order for like charged and all charged particles are given, 
respectively, by 

{n+{n+ - 1)) 



= 2 (for like sign particles) , (31) 



and 



K 



((n+ + n_)^ - (n+ + n_)) 



(n+(n+ - 1)) + 
2(n+)2 

= 1.5 (for charged particles) . (32) 
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For y^=34.5GeV, one can investigate the fluctuations in a single cell, i.e. fluctuations 
induced by pure Bose-Einstein distributions(cf. eq ( ^2|) ) at M=5 since Sy=0.8 and 
AF=4.0. Our simulation gives Fj''(M = 5) ^ 1.27, and Fj"(M > 16) ^ 1.38. The 
reduction from 1.5 can be attributed to the smearing due to varying cell location(cf. the 
item 1 of the step 2 of sec. 11 A ) and also to the global charge conservation. In fact, we 
obtain F^^'^M = 5) :^1.32, and F'^{M > 16) ~ 1.46 when all the conservation laws are 
not imposed. 

Now the results are compared with experimental data|jl6|, |2l| for both exclusive and 
inclusive factorial moments in Figs. 9 and 10. Our results approximately reproduce the 
experimental data. In particular, the agreement is very good for ^/s = 34.5 GeV. In our 
model, the parameter 6y plays a crucial role in reproducing the bin size dependence of the 
factorial moments. As M increases from unity, the factorial moments increase due to the 
BEG and then start to saturate for M > AY/ 6y^. In this sense, the parameter Sy in our 
model may be interpreted as an effective correlation length in the rapidity space due to the 
BEG. Our results show that the BEG considerably contribute to the increase of the scaled 
factorial moments(the so-called intermittency) . However, at y/s = 91.2 GeV our results 
of both exclusive and inclusive factorial moment systematically underestimate compared 
with experimental data for all orders. This underestimation becomes remarkable as the 
order of the moment increases. It may be due to three jets events at high energy not 
taken into account in our calculation. At LEP energy, events with three jets take place 
with a signiflcant probability. 



Figure 9. (a)-(c) 
Figure 10. (a)-(c) 



4.5 Bose-Einstein correlations 

Here the results for the Bose-Einstein correlation functions are compared with data at 



y/s = 34.50, and 91.2 GeV [0, |22[. The following correlation functions C{Q) (or 
C{Q^)) are used for this analysis: 

CiQ) = Ar±±(Q)/iV-/(Q) or CiQ^) = iV±±(g2)/Ar-^(Q2), (33) 

where N^^{Q) (or N^'^{Q'^)) is the number of like sign pairs with four momentum dif- 
ference: = {pi — OT Q = \fQ^. The so-called reference sample, the denominator 
N'^'^^ is the number of pairs without Bose-Einstein correlations. TASSO and ALEPH 
collaboration use the number of unlike sign pairs as N"^^^ . 'Mixed pair reference' is also 
used by ALEPH. In this case, a pair is taken from different events. Our results of both 
the TT- and the p- models are shown in Figs. 11 and 12, respectively. We also present 
energy dependence of the BEG correlation function C{ff) in Figs. 13(a) for the 7r-model 
and (b) for the p-model. 

TASSO Collaboration compared their data[^ (exclusive factorial moments at ^/s = 34.5 GeV ) 
with results of various event simulators of multiparticle production model: Webber, Hoyer and Lund 
model version 6.2 and 6.3 (see Refs in Those models can also roughly reproduce the experimental 

data. The BEC are not taken into account while cascade processes are included in them. 
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Figure 11. (a),(b) 
Figure 12. 



Our model produces a characteristic enhancement of C{Q) {C{Q'^)) in small Q (Q^) 
region. The result of the p-model is in better agreement with experimental data than the 
TT-model for the whole region of Q^. Wc found that the enhancement of the correlation 
function C{Q) (C(Q^)) strongly depends on the size of Sy. When Sy becomes large, the 
slope of the correlation function becomes small. If Sy is larger than the rapidity interval 
for the full phase space, i.e. 6y > 2ymax, the correlation function has no Q dependence, 
C{Q) = 2. On the contrary when 6y approaches to zero, the slope of the correlation 
function becomes large and C{Q) approaches to unity for Q > 0. 

Although our model has no explicit information on the space-time structure of particle 
emission points, we can extract the size and the lifetime of the particle sources by fitting 
a theoretical formula to the result with mixed pair reference from our simulation. Here 
we use the following fitting formula for Q >5 MeV/c: 

CiQ) = c[l + Acxp(-i?g)], (34) 
C{Q^) = c[l + \exp{-R^Q% (35) 

where c is the normalization constant, A is the so-called chaoticity parameter, and R is the 
size parameter. Wc don't use data for Q <5 MeV/c in the fitting because the statistical 
errors are very large there. 



Figure 13. (a) and (b). 



In Table El, the results of the fitting at various energies arc summarized for both the 
TT-model and the p-model. With increasing energy of the system, R decreases while A 
tends to increase. We find that R in the p-model is larger than that in the TT-model. We 
also try to extract the longitudinal size Rl, transverse size Rt, and lifetime r by using 
the following formulae in our fit:. 

C{Ql) = cl[1 + ALexp(-i?iQi)], C{Qt) = ct[1 + ATexp(-i?^g|)], 

C(A£;) = c^[l + A^exp(-T^A£;2)], (36) 

where Q\ — \pl\ — PL2p, Qt — \Pti — Pr2p and AE = |ei — €2]. The results are presented 
in Table IV. The source size decreases with increasing ^/s, while the life time is almost 
independent of ^/s. It is reasonable that C{Qt) has no strong (^T-dependence because 
the BEG is not taken into account in the transverse momentum space in our model. 



Table HI and IV 



The value of the normalization constant c may be understood as follows. The num- 
ber of unlike charged pairs and that of (++) pairs for large Q may be proportional to 
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(statistical) combinatorial numbers (n_|_n_) and {n+{n^ respectively. Therefore the 

asymptotic value of the correlation function may be estimated as 

lim C{Q) = c (37) 
_ (n+(n+ - 1)) 



(St 



where we have used charge conservation in deriving the second line. Eq. (^Hj) with the ap- 
proximation (n^) ~ {n+)^ and the experimental values of {nch){= 2(n+)) yields c=0.785, 
0.855 and 0.901 for ^5=14.0, 34.5 and 91.2 GeV,respectively in qualitative agreement 
with the result given in Table IE. 

We also investigate the influence of the charge and energy-momentum conservation on 
the Bose-Einstein correlations. Both energy-momentum conservation and charge conser- 
vation reduce the values of the correlation functions in the whole Q (or Q"^) range as shown 
in Fig. 14. The effect of charge conservation is stronger than that of energy-momentum 
conservation. If "no conservation" events are used, C{Q) approaches to a larger constant 
as Q increases. 



Figure 14. 



As mentioned in Section 111 B, the BEG in the p-model is imposed for like-charged 
pions by using the enhancement factor. By changing this enhancement factor, we can 
control the strength of the BEG. For example, one can reduce the BEG substantially by 
putting nf = = in eq.(|25|) as shown in Fig. 15. 



Figure 15. 



4.6 2-particle rapidity correlations 

We investigate the 2-particle correlations in terms of rapidity variables in detail. The 
experimental data reported by TASSO Gollaborationfl^ at ^3=34. 5 GeV is used in 
comparison with our calculations. When the 2-particles of type "a" and "6" ( for example, 
a,b = TT"*", 7i~) have rapidities yi and y2, respectively, the two particle rapidity correlation 
function R'^''^\yi,y2) is defined as follows: 

'^"^^> = /P«(yOp<"(y.)"'' 

where p^"''''\yi,y2) and p^"'\yi) is the two particle density and the one particle density, 
respectively: 

P^''''\yi,y2) = (40) 
p'^'Hyi) = • (41) 

a ay 
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Here a is the total cross section and / is the normahzation constant given by 



< na{nb - Sab) > 
<na> <nb>' 



(42) 



Jab 



where 

for a=b 

for aT^b (43) 
when particle type is not distinguished, 

and Ua and rib are the multiplicities of "a" and "6" ,respectively. 

Instead of using R2''^\yi,y2), in practical measurements by TASSO Collaboration, 
they use 



r 1 






R2'''\yi,ytrigg) = / ^ dyt i^^^'^^^/i, ^t) , 

y/3 — ya Jy<x 



(44) 



where ytrigg refers to a rapidity interval {ya, yp)- These trigger rapidity intervals are trig- 
ger 1 ( -5.50< y < -2.50), trigger E ( -2.50< y < -1.50), trigger H ( -1.50< y < -0.75) 
and trigger IV ( -0.75< y < -0.00). First we investigate a case where sign of charge is 
not distinguished and hence 5ab—^- The results are compared with experimental data in 
Fig.l6(a)-(d). The experimental data are approximately reproduced by the simulation for 
every trigger regions. Although our model predicts a negative correlation for large y {y >2 
) in the trigger I, a positive correlation is reported by TASSO. For the trigger regions E, 11 
and IV, calculated correlation functions show positive bumps around the trigger rapidity 
interval in agreement with experiment. The effect is apparently caused by the BEC. 



Figure 16(a)-(d) 



For confirmation, we calculate the 2-particle rapidity correlation R^^ {y , ytrigg) and 
{y : ytrigg) ■ See Fig.l7(a)-(d). The clear positive correlations are seen near the trigger 
rapidities in the correlation function R^^^\yi,y2). On the other hand, such positive 
correlations are absent in i?*^^ \yi,y2)- From these investigations, one can conclude 
that the BEC (identical particle effect) plays a crucial role in an enhancement of the 
R^^{y, ytrigg) near the trigger rapidity ytrigg- The correlation functions are also affected 
by the effects of p meson production and their decay. Results of the p-model are shown in 
Fig. 18. It is found that the p-model gives better agreement with experimental data than 
the TT-model for all trigger rapidity intervals. 



Figure 17 (a)-(d) 



There are some discrepancies between the model result for R{y, ytrigg) and the data 
as most evidently seen in Fig. 16(a). The situation is qualitatively the same in both the 
TT- and the p- models. This is an indication that there are unknown long range as well as 
short range correlations besides BEC, e.g. those due to higher resonances and/or three 
jets. 
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Figure 18. (a)-(d) 



In analogy to C {Q) , another 2-particle rapidity correlation function can be defined by 
using the rapidity deference Ay = — 2/2! in place of Q: 

C{Ay) = Z , . • (45) 

This function allows us to measure the size of the model parameter 5y. We expect that 
slope of C{Ay) will change at Ay ^ 5y. In this paper, we assume that 5y is independent 
of y. This assumption can be checked by measuring C{Ay) in various rapidity regions. 
We strongly urge experimentalists to measure C{Ay). Our prediction for this correlation 
function is shown in Fig.l9(a)-(c). 



Figure 19. 



4.7 Forward-backward correlations 

Finally we compare the results of the forward-backward correlation with experimental 
data|TE|, |2^. It is known experimentally that there are weak (positive) correlations be- 
tween the multiplicity in the backward hemisphere and the mean multiplicity (np) 
in the forward hemisphere. The direction of the "forward" or "backward" is assigned at 
random in the simulation. The comparison of our results with experimental data is shown 
in Fig.20. It is found that our model reproduces well experimental data for every ^/s. 
This correlation is usually parameterized by a linear form: 

{np) = a + b ub- (46) 

The values of the parameters a and b extracted from the vr-model and the corresponding 
experimental values are given in Table V. 



Figure 20. 
Table V 



The correlation strength b depends on the size of 6y. When Sy becomes large, the 
correlation strength b also becomes large. When the multiplicity distribution is extremely 
narrow, the correlation strength b is negative. For example, if the multiplicity distribution 
is delta function type, P{nch) = 5{nch — {nch)), the correlation strength b is —1. Note that 
the dispersion of the multiplicity distribution is also sensitive to the size of 6y. Increase 
of Sy enhances the dispersion of the multiplicity distribution. This enhancement of the 
dispersion also increases the correlation strength b. We found that the forward-backward 
correlation is not significantly affected by the energy-momentum conservation. 
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5 Conclusions and discussions 



We have investigated various correlations and fluctuations observed in multiparticle pro- 
duction in e^e~ annihilation by using a new statistical model constructed on the basis 
of the maximum entropy method. This model allows us to investigate the crucial roles 
of the Bose-Einstein Correlations (BEC) in multiparticle production phenomena. The 
model has three parameters T, and 5y. The "partition" temperature T and the "chem- 
ical" potential ji characterize the single particle spectra. The BEC is characterized by 
the third parameter 8y. Based on this model we have constructed two kinds of event 
generators which correspond to two extreme cases. One case is the 7r-model which as- 
sumes that all pions are produced directly, i.e., e"'"e~ — > tt tt • • • . The other is the 
p- model where it is assumed that all pions are produced via the decay of p mesons, i.e., 
e'^e~ p p > TTTT In the course of the event generation, energy-momentum 
and charge conservations are imposed by event selection. Thus we can study the effects 
on single particle spectra and correlations caused by these conservation laws. We found 
that some observables are significantly affected by them. 

Various kinds of the experimental data on single particle spectra and many particle 
correlations are well reproduced by our model. We have compared the results obtained 
from our simulations with experimental data on multiplicity distributions, rapidity distri- 
butions and rapidity dependence of the mean transverse momentum. For many particle 
correlations, we compare our calculations with experimental data on factorial moments, 
Bose-Einstein correlations, 2-particle rapidity correlations and forward-backward corre- 
lations. We have also calculated the charge correlations in the rapidity space. We hope 
experimentalists will measure this observable. 

Our results for all observables on correlations, including the dispersion of the multi- 
plicity distributions, exhibit a strong 5y dependence. We found that 5y plays an essential 
role in explaining the behavior of some 2-particle correlation functions, scaled factorial 
moments and so on. Once the value of the parameter 5y is determined by fitting our result 
on the rapidity dependence of the dispersion of the multiplicity distribution to that of the 
experimental data, other correlation and fluctuation data are systematically reproduced 
by using the same 5y value. The parameter 5y may correspond to an effective correlation 
length (in the rapidity space) of the BEC. Constancy of the size of the 5y in the whole 
rapidity space is the most characteristic point in our model. This simple assumption may 
be checked by measuring the correlation C(A|/) at various rapidity regions. There are 
possibilities that it depends on the rapidity y and/or the multiplicity. 

In Bose-Einstein correlations and 2-particle rapidity correlations, in particular, we 
observe a clear difference between the prediction from the 7r-model and that from the p- 
model. We found that the p-model gives better agreement with those experimental data 
than the 7r-model. This means that two-particle correlations are sensitive to production 
and decay processes of resonances and experimental data on those observables suggest 
that there are significant contribution from resonance decay. No significant difference 
between the vr-model and the p-model is observed in the single particle spectra and other 
correlation data. 

Finally, we would like to point out two important aspects of our model. First, we 
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would like to emphasize that the BEC are incorporated on event-by-event basis in our 
model. It thus provides a useful theoretical tool for event-by-event analysis of the BEC. 
This is an important feature of our model not shared by any other event generators. 
Second, information on the space-time structure of multiparticle production is apparently 
not contained in our statistical model. Nevertheless, it gives the BEC from which one 
can extract the information on the "apparent" source size and the chaoticity. It thus 
appears that the information on the space-time structure is "hidden" in our model. The 
fundamental parameter 5y may have some relevance to the question. 
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Table captions 

Table I: Best fit values of /i, T and 5y for various with the observed values of {rich) 
and (pt) used in the 7r-model. 

Table II: Values of the parameters in the p-model. 

Table HI: The results of i?, A and c in the tt- and the p- models. 

Table IV: The results of Rl, Rt and r in the 7r-model. Corresponding results of the 
chaoticity parameters Al, At and A,- are also presented. 

Table V: The results of the parameters which characterize the forward-backward cor- 
relation. The experimental values are given in the parentheses. 

Figure captions 

Fig.l: Total energy of events generated for the mean total energy (E^tot) = -\A = 34.5 
GeV and the energy window of the full width 25E=QA^. 

Fig. 2: Comparison of the results of the 7r-model with the experimental data on the ycut 
dependence of DlJ {rich) for (a) ^^U.Q GeV, (b) 34.5GeV and (c) 91.2 GeV. 

Fig. 3: Comparison of multiphcity distributions in the 7r-model for (a) ^/s—lA.{) GeV, 
(b) 34.5GeV and (c) 91.2 GeV. 

Fig. 4 (a): KNO scaling of the multiplicity distributions in the tt- model at energies 
-y/s=14.0, 34.5 and 91.2 GeV. (b): Comparison of the tt- model(sohd line) with the 
experimental data for energy GeV with ycut ~ 6.5(full phase space), 2.0 

and 1.0. 
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Fig. 5: Effects of conservation laws on ycut dependence of D^j^/ (rich)^- 

Fig.6 (a): Prediction of charge correlation DlJ2Dl for a/s = 14.0, 34.5 and 91.2 GeV. 
(b): Effects of conservation laws on charge correlation L>^^/2D^. 

Fig.7: Rapidity distribution dN/dy for = 14.0, 34.5 and 91.2 GeV in the tt- 
model(open squares). Experimental data for y/s — 14.0 and 34.5 GeV are shown 
by solid circles. 

Fig.8 (a): Rapidity dependence of the mean transverse momentum (pt) for ^/s — 14.0, 
34.5 and 91.2 GeV. (b): Effect of the thrust axis and/or the energy-momentum 
conservation for energy ^/s — 34.5 GeV. 

Fig.9: The exclusive factorial moment for (a) ^/s — 14.0, (b) 34.5 GeV and (c) 91.2 
GeV in the 7r-mo del (solid line). Experimental data are shown by open circles. 

Fig. 10: The same as in Fig.9(a)-(c) for the inclusive factorial moments in the 7r-model. 

Fig. 11: Comparisons of the results for C((5^) from the tt- and p- model with the 
experimental data for (a) ^/s — 34.5 GeV and (b) ^/s — 91.2 GeV. The number of 
unlike pairs is used as the reference. 

Fig. 12: The same as in Fig.ll for ^/s = 91.2 GeV with the mixed pairs as the reference 
sample. 

Fig. 13: Energy dependence of C{Q) with mixed pair reference in the (a) tt- and (b) p- 
models. 

Fig.l4: Effect of conservation laws on the BEG correlation function C{Q). 
Fig. 15: The role of the enhancement factor in the p- model. 

Fig. 16 (a): Comparison of the result from the 7r-modcl with experimental data on 
the 2-particle rapidity correlation function R{y, ytrigg) at ^/s— 34.5 GeV for (a) the 
trigger I [-5.50,-2.50], (b) the trigger I [-2.50,-1.50], (c) the trigger HI [-1.50,-0.75] and 
(d) the trigger IV [-0.75,0.00]. 

Fig. 17: The 2-particle rapidity correlation functions R^'^iy, ytrigg) and R^^iy, ytrigg) 
in the 7r-model for ^/s— 34.5 GeV with (a) the trigger I [-5. 50,-2. 50], (b) the trigger 
I [-2.50,-1.50], (c) the trigger HI [-1.50,-0.75] and (d) the trigger IV [-0.75,0.00]. 

Fig. 18: Comparison of the results of p- model with experimental data on the 2-particle 
rapidity correlation function ytrigg) at \/s— 34.5 GeV for (a) the trigger I [- 
5. 50,-2.50], (b) the trigger I [-2.50,-1.50], (c) the trigger HI [-1.50,-0.75] and (d) the 
trigger IV [-0.75,0.00]. 

Fig.19: Predictions for the 2-particle rapidity correlation function C{Ay) for ^/s—14:.0, 
34.5 and 91.2 GeV. 

Fig. 20: Comparison of the forward-backward correlation in the 7r-model with experi- 
mental data for ^/s^U.O, 34,5 and 91.2 GeV. 
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[GeV] 


5y 


/i [GeY] 


T [GeV] 


{rich) 


{Pt) [GeV] 


14.0 


0.60 


-1.71 


2.47 


9.30 


0.334 


22.0 


0.70 


-2.43 


4.26 


11.30 


0.377 


34.5 


0.80 


-3.08 


5.75 


13.59 


0.422 


43.6 


0.85 


-3.77 


7.83 


15.08 


0.446 


91.2 


1.20 


-4.73 


15.4 


20.80 


0.521 



Table I 



V~s [GeV] 


hp 


hn 


A^P [GeV] 


Tp [GeV] 


14.0 


1.20 


0.60 


-3.61 


3.05 


34.5 


1.60 


0.80 


-10.64 


9.00 


91.2 


2.40 


1.20 


-17.34 


20.27 



Table H 



[GeV] 


model 


function 


i?[fm] 


A 


c 


14.0 


71 {Sy=0.6) 


Gaussian 
Exponential 


0.56±0.02 
0.72±0.04 


0.46±0.02 
0.66±0.02 


0.84±0.00 
0.83±0.00 


, (%=l-2) 
P {Sy„=0.6) 


Gaussian 
Exponential 


0.77±.04 
1.17±0.09 


0.28±.01 
0.43±0.03 


0.85±.00 
0.85±.00 


34.5 


n iSy=0.8) 


Gaussian 
Exponential 


0.50±0.02 
0.62±0.02 


0.58±0.02 
0.83±0.02 


0.90±0.00 
0.88±0.00 


(H=l-6) 


Gaussian 
Exponential 


0.61±0.02 
0.80±0.04 


0.39±0.01 
0.54±0.02 


0.91±0.00 
0.91±0.00 


91.2 


TT (5y=1.2) 


Gaussian 
Exponential 


0.36±0.01 
0.42±0.01 


0.62±0.02 
0.92±0.01 


0.96±0.01 
0.91±0.00 


, (%=2.4) 


Gaussian 
Exponential 


0.48±0.01 
0.63±0.02 


0.34±0.01 
0.49±0.01 


0.96±0.00 
0.95±0.00 



Table HI 
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V^s [GcV] 


RlM 


AL[fm] 




AT[fm] 


r[fm] 


A. 


14.0 
91.2 


0.55±0.03 
0.34±0.01 


0.32±0.02 
0.34±0.01 


0.29±0.16 
0.13±0.05 


0.07±0.03 
0.04±0.01 


0.20±0.00 
0.20±0.01 


0.20±0.01 
0.13±0.00 



Table IV 



x/i [GeV] 


a 


b 


14.0 ((5y=0.60) 


4.36 


0.054 (0.085± 0.014) 


31.5 (()y=0.80) 


6.18 


0.071 (0.089± 0.003) 


91.2 (5y=1.20) 


9.03 


0.137 (0.118± 0.009) 



Table V 
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